Formation and decay of electron-hole droplets in diamond 
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We study the formation and decay of electron-hole droplets in diamond at low and high tem- 
peratures under different excitations by master equations. The calculation reveals that at low 
temperature the kinetics of the system is similar to that of a direct semiconductor, whereas at high 
temperature it is metastable and similar to an indirect semiconductor. Our results at low temper- 
ature are consistent with the experimental findings reported by Nagai et al. [Phys. Rev. B 68, 
081202 (R) (2003)]. The kinetics of the e-h system in diamonds at high temperature under both 
low and high excitations is also predicted. 



PACS numbers: 71.35.Ee, 71.35.-y, 81.05.Uw 

Photoexcited electron-hole (e-h) systems in semicon- 
ductors provide a unique opportunity to understand 
quantum many-body phenomena with Coulomb interac- 
tions. In the dilute density region, an electron and a hole 
are combined to form a neutral bound state — an exciton. 
At low temperature, a dense exciton gas condenses into a 
Uquid phase with a metaUic character in the form of e-h 
droplets (EHDs). This macroscopic metallic phase has 
been extensively investigated in the past three decades 
in indirect semiconductors such as Ge and Si.^ The tran- 
sition between EHD and exciton gas is considered anal- 
ogous to a classical liquid transition in water, and the 
EHD formation is well understood with a classical nu- 
cleation theory.^ However, the formation and decay of 
EHD in photoexcited semiconductors are not only deter- 
mined by the collection and evaporation of excitons on 
the surface of EHD but also by carrier recombination. 
In direct semiconductors in particular, a fast recombina- 
tion process overcomes the thermal kinetics of carriers. 
Therefore e-h pairs annihilate before a small e-h ensemble 
grows to become a macroscopic-size EHD. Consequently, 
the phase transition is shown to be of second order. 

This competition between thermal kinetics and recom- 
bination of carriers is also apparent in indirect semicon- 
ductors. In traditional indirect semiconductors such as 
Ge and Si, at a certain high temperature the evaporation 
rate is larger than the recombination rate, which makes 
the kinetics of EHD formation similar to that of classical 
nucleation. In this case, EHD formation exhibits a hys- 
teresis effect and the average drop size is large. However 
when the temperature is sufficiently low, the thermal ki- 
netics is suppressed. The dominant recombination effect 



makes the e-h system behave like those in direct semi- 
conductors, i.e., no hysteresis effect and a small average 
number of pairs per cluster (ANPC). Under this condi- 
tion the exciton-EHD phase transition changes from first- 
to second-order.'^''* This density and temperature region, 
where the thermo-dynamical phase diagram is distorted, 
is attractive to scientists because the quantum statistics 
of the quasi-particles is dominant, and a hidden collec- 
tive phase including Bose-Einstein condensation might 
appear. In order to understand such a rich variety of 
macroscopic phases, it is important first to evaluate the 
kinetics of the liquid-gas transition in the photoexcited 
indirect semiconductors in the low temperature region. 
Nevertheless it is not very realistic for conventional in- 
direct semiconductors because of their narrowness in en- 
ergy scale. 

Diamond is a wide band gap indirect semiconductor 
with a band structure similar to those of Ge and Si, and 
is a good candidate to study carrier dynamics. More- 
over, because of the small dielectric constant of diamond 
the screening of the Coulomb interaction between carri- 
ers is small. Thus one can treat e-h system in diamond 
in wide energy scale. Recently Shimano et al. evalu- 
ated the character of EHD in diamond by time-resolved 
luminescence measurements and reported a higher criti- 
cal temperature, larger work function, larger density, and 
shorter lifetime for EHD in diamond compared to Ge and 
Si.^ Consistent values are also obtained from an analysis 
of the luminescence spectra under quasi-cw excitation.^''' 
The dynamics of the EHD formation at 12 K under a dif- 
ferent excitation density has also been studied by Nagai 
et al.^ It was observed that after photoexcited carriers are 
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cooled rapidly into a supersaturated exciton gas within 
several tens of picosecond, spatial condensation of dense 
exciton gas into EHD occurs within a few hundred pi- 
coseconds. In this report we investigate theoretically the 
kinetics of EHD formation and decay in diamond. First 
we use a discrete master equation theory developed by 
Haug and Abraham^° to investigate the femtosecond ex- 
citation in diamond at a low temperature regime where 
only small c-h clusters are formed. We then use the con- 
tinuous master equation theory developed by Silver^^ and 
by Koch and Haug^^ to investigate the dynamics at high 
temperature regime where the average drop size is too 
large to be treated discretely. Finally we compare our 
results with the experimental measurements by Nagai et 
al.^ The division between the low and high temperature 
regimes in diamond is ~60 K where the thermal evapo- 
ration rate equals to the recombination rate.^"^ 

For a discrete master equation formalism, if the con- 
centration of clusters containing n e-h pairs at time t is 
denoted by f{n,t), the master equation describing the 
evolution of f(n, t) is: 

d 

-g^f{n,t)=jn-i-jn (1) 

for n > 2, where j„ is the net probability current between 
the clusters with n and n+1 e-h pairs: 

jn = 9nf{n, t) - ln+if{n +l,t) . (2) 

In this equation l„ and <?„ are the gain and loss rates of a 



cluster with n e-h pairs. The gain rate is obtained from 
the assumption that excitons with a density are col- 
lected at the surface of a cluster, and is approximated by 
9n = bfixfi^/^ with b = AttRIvx- Rq = (3/47rpo)"^/^ is the 
Wigner-Seitz radius of the EHD and — y/ kT/2TTmx 
is the thermal velocity of excitons with an effective mass 
rux- po denotes the EHD density. The loss rate is com- 
posed of the sum of the evaporation rate a„ and the re- 
combination rate n/r„. In = an + n/Tn- The evaporation 
rate is given by a time-independent Richardson-Dushman 
current, a„ = bDx exp[{-(j> + can~'^/^) / {kT)]n'^/^ , where 

the effective density of state 
of exciton, jx is the degeneracy of the exciton ground 
state and ccrn"^/^ represents the correction of the bind- 
ing energy due to surface effect with a denoting the sur- 
face energy of the EHD. These equations are solved to- 
gether with the continuity equation: 

l.m,t) = Git)-f:^-2n-±j., (3) 

71=1 " n=2 

with G{t) representing the excitation pulse. 

G{t) = Goe-(*-*°)'/*? (4) 

with tp standing for the width of laser pulse which is 0.1 

ps throughout this paper. It is noted that the coalescence 
of clusters larger than excitons is neglected. 



TABLE I: EHD and exciton parameters for diamond which arc used in the calculation 





symbol 


value 


unit 


Ref. 


mean EHD lifetime 


Td 


1 


ns 


5 


mean exciton lifetime 




100 


ns 




work function of EHD 




50 


moV 


5,7,15 


surface energy of EHD 




1.2 


erg/cm^ 


16 


a(T) = ao(l - (T/Te)^) 


Tc 


165 


K 


5 


exciton degeneracy 


7 


12 






effective mass of exciton 


rrix 


7.92 


10-^1 kg 




e-h density of EHD 


Po 


1.0 


10^° cm-3 


5 



For the case of high temperature where n is too large 
to be treated discretely, we turn to the equation of 
moments. The z^-th moment of EHD distribution is de- 
fined as 

oo 

x^{t) = j rff{n,t)dn (5) 

Tie 

in which is a critical size where the stationary dis- 
tribution has a minimum. All clusters smaller than ric 
are counted as excitons, while clusters larger than rig are 



treated as EHDs. Under this approximation, the "exci- 
ton density" is given by: 

nx{t) = jnf{n,t)dn. (6) 
1 

Here Uc is calculated approximately by equating loss rate 
and the gain rate with the recombination loss neglected: 

^ nx{t)bny\ (7) 



3 



The equations of moments are given by: 



(8) 



d d X 

^•^ = '^c^^O - + X^_i/3b{ns - Ua;)] , (9) 

at Td 



dt 



with Td denoting the mean EHD Hfetime and Ug be- 
ing the saturated exciton density, Ug = Dxexp{—(j)/kT). 
V = 1/3, 2/3, 1, 4/3, 5/3, and 2. Finally the continuity 
equation is 



dt 



Tl Qj 

~ - Uc-^xo + X2/zh{ns - Tlx) + G{t) , (10) 



dt 



in which reprcscints the mean exciton lifetime. The 
second term describes the change in EHD density and 
the third term is due to free exciton evaporation and 
collection by EHD. Equations (8) to (10) form a closed 
set of equations. The expressions of and f{nc,t) in 
Eq. (8) are given in Appendix A. 

We first study the kinetics of the c-h system in diamond 
at 12 K by using the discrete master equations, under the 
experimental conditions similar to the ones used by Nagai 
et al.^ The material parameters are listed in TABLE 1.^^ 
The results of our calculation are plotted in Figs. 1-4. 

In Fig. 1 we present the kinetics under an excita- 
tion with Go = 6.6 X 10^^ cm~^s~^ and tp = 0.1 ps 
which corresponds to a ~0.2 mJ/cm^ excitation in the 
experiment. ^■"'^^ The cluster concentration versus time 
and the number of e-h pairs in a cluster as well as 
the concentration of some selected clusters versus time 
are plotted in Figs. 1(a) and (b) respectively. The fig- 
ures indicate that the concentration of small clusters rise 
faster than that of large cluster because the current j„ 
in Eq. (1) flows from small clusters to larger ones. It is 
also seen from the figures that the system reaches quasi- 
equilibrium at about 120 ps and this quasi-equilibrium 
lasts about 200 ps. 

In the experiment the peak energy of the EHD emis- 
sion band shifts toward the low energy side during first 
200 ps,* which suggests that large clusters are formed at 
a longer time. Meanwhile from the fact that there is little 
change of the luminescence in the low energy regime, we 
conclude that the rate of the formation of small clusters 
from excitons and the rate of coalescence to large clus- 
ters are nearly the same. As a result the concentration 
of small clusters rises and quickly reaches a steady value. 
These features are consistent with our calculation. The 
luminescence indicates that the system reaches its quasi- 
equilibrium in about 200 ps, which is comparable to the 
120 ps value we obtained. It is also seen in Fig. 1(b) that 
excitons (n = 1) decay much slower than clusters with 
n > 2 — because the lifetime of excitons Tx is larger than 
that of e-h pairs in e-h clusters and also because exci- 
tons get additional compensation from the recombination 
of large clusters. 

In order to monitor the average size of the cluster, we 
plot the ANPC (n) = Er=i 0/ Er=i /K as a 
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FIG. 1: (a): Cluster concentration versus time and number of 
e-h pairs per cluster for diamond at T = 12 K under Gaussian 
pulse excitation with Go = 6.6 x cm~^s~^ ,tp = 0.1 ps, 
and to = 1.0 ps. (b): Concentration of some selected clusters 
(solid curve) and ANPC (n) (dashed curve) versus time. 



function of time in Fig. 1(b). The ANPC is less than 
3, and is different to that obtained for direct band-gap 
semiconductors. ^°'^* Note that in the present model the 
coalescence of clusters larger than excitons {e.g. biexci- 
tons) is neglected. The coalescence of clusters adds to 
the cluster formation mechanism, thus this approxima- 
tion leads to a smaller ANPC in our results. Despite 
this, the ANPC is still too small to residt in the forma- 
tion of macroscopic EHDs. The smallness of ANPC and 
the shortness of the time during which the system is in 
quasi-equilibrium indicate that the system is characteris- 
tic of a nonequilibrium, similar to an e-h system in direct 
semiconductors . 

We also show the time evolution of the exciton den- 
sity nx and the total density of all e-h pairs condensed in 
clusters larger than excitons: = Ej^2 ^/('^j ^) Fig. 
2. Thus we can compare the evolution of these densities 
with that of excitons and the integrated EHD lumines- 
cence intensities. We can see that excitons slowly con- 
dense into EHDs and na reaches a maximum around 150 
ps, which corresponds to the 260 ps experimental data. 

There arc some differences between these densities and 
the luminescence intensities in the experiment. First, the 
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FIG. 2: Exciton density Ux (solid curve) and the total den- 
sity of e-h pairs in all clusters larger than exciton rid (dashed 
curve) for diamond at T = 12 K under the same excitation as 
in FIG. 1. 



times when exciton and EHD densities reach their max- 
ima are about one half those in the experiment. This 
originates partly from the simplified excitation model we 
use. In this model the only excitations generated by 
the laser pulse arc excitons. In reality, these excitations 
should be e-h plasma and they are always overheated. Af- 
ter the excitation the e-h system is cooled down in several 
tens of picoseconds.^ This relaxation process affects the 
kinetics of exciton. The absence of this process results in 
a shorter formation time in the calculation. Second, at 
the equilibrium stage 90% of the excitation is converted 
into the EHD phase, while in the experiment only about 
50% is converted. This may be partly due to the diffusion 
of the e-h pairs, and therefore the actual exciton density 
is smaller than our evaluation. Nevertheless the diffusion 
effect might be marginal since the time scale is only 300 
ps and the initial carrier distribution with a penetration 
depth of 15 /zm is less spatially inhomogeneous. It is re- 
ported that the effect of diffusion of an e-h system in Si is 
negligible during the 200 ps after excitation.^'' Moreover, 
the high density of e-h pairs causes large Auger recombi- 
nation in EHD and a repulsion of excitons from EHD by 
the phonon wind.^'^ Thus the efficiency of the collection 
of excitons which collide with EHD decreases, and less 
excitons are converted into EHD. These effects also slow 
down the formation process as the gain rate gf„ in Eq. 
(2) is proportional to exciton density, and thus lead to 
a longer formation time in the experiment compared to 
our results. 

We now discuss kinetics under a higher excitation. Fig. 
3 shows the same calculation as in Fig. 1 but with a much 
higher excitation, i.e., ninety times as large as in Fig. 
1. This intensity corresponds to an excitation of ~17 
mJ/cm^ in the experiment.^ It is seen from the figure 
that compared to the case of low excitation in Fig. 1 , the 
exciton density decays during the first 20 ps. This ex- 
plains the absence of exciton luminescence in the exper- 
iment performed at a similar excitation.^ It is also seen 
that although the peak exciton density is almost pro- 
portional to the excitation intensity, the exciton density 
in the quasi-equilibrium is smaller for larger excitation 
intensities. In Fig. 1(b) and Fig. 3(b), the exciton den- 
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FIG. 3: (a): Cluster concentration versus time and number of 
e-h pairs per cluster for diamond at T = 12 K under Gaussian 



pulse excitation with Go = 5.6 x 10 cm s 



in 



0.1 ps, 



and to = 1.0 ps. (b): Concentration of some selected clusters 
(solid curve) and ANPC (n) (dashed curve) versus time. 



sity in quasi-equilibrium is about 10 cm" under an 
excitation of 0.2 mJ/cm^ and only 10^^ cm"^ under 17 
mJ/cm^. Moreover the ANPC increases very little de- 
spite such a large increase in excitation intensity. This 
behavior confirms the second-order nature of the exciton- 
EHD transition at this low temperature, and it is similar 
to what was discovered in Ge and Si at sufficiently low 
temperatures.'! j^- yndej-stood that the larger excita- 
tion intensity makes the concentration of clusters grow 
up more rapidly since the gain Qn in Eq. (2) is propor- 
tional to the exciton density. 

Before we discuss the high temperature case, we an- 
alyze the dependence on the excitation intensity at low 
temperature. In Fig. 4 we compare the time delays at 
which Ud reaches the maximum and half maximum as 
well as the time delay needed for the exciton density 
to reach a maximum under different excitation intensi- 
ties [Fig. 4(a)] with those in the experiment [Fig. 4(b)]. ^ 
The figure shows that our results are in good qualita- 
tive agreement with experiment. In both figures the for- 
mation time of e-h cluster decreases with the increase 
of the excitation intensity, and the time when the ex- 
citon density (luminescence intensity) reaches its max- 
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FIG. 4: Calculated (a) and experimental (b) e-h droplet for- 
mation times and exciton formation times versus the excita- 
tion intensity. •: time used when rid reaches its maximum; 

time used when Ud reaches its half maximum; A: time 
used when Hx reaches its maximum. 



imum is independent of the excitation intensity. One 
may find that the time of the maximum density of ex- 
citon in our results is nearly zero which is smaller than 
the experimental value around 20 ps. This is due to the 
fact that we use a simple model of Gaussian excitation 
in which the relaxation process of the e-h system is ne- 
glected. As mentioned before, the relaxation, diffusion, 
and the Auger processes slow down the EHD formation 
process. Thus our calculations produce a relatively small 
formation time. 

We now study the kinetics of diamond at high temper- 
ature (100 K) using equations of moment Eqs. (8-10). It 
is noted that this method only allows one to study the 
formation and decay processes close to the steady state. 
We follow the path used for the low temperature case, 
i.e., first a small excitation and then a high excitation. 

In Fig. 5(a) we present the evolution of the "exciton 
density" Eq. (6), the three integer moments of the EHD 
distribution Eq. (5), the average size of EHD n = xi/xq 
and the relative mean square of EHD distribution 



(An)2 = [X2/X0 - (xi/xo) ]/{xi/xo) 



(11) 



under an excitation of 1.8 mJ/cm^. Very different from 
the results at low temperature where only small clusters 
are formed, here one finds that the average drop size is 
very large: about 10^. Nevertheless the density of all e-h 
pairs which are condensed in EHDs, i.e. xi, is rather 
small, less than 10^^ cm~'^. This can be understood eas- 
ily: For a high temperature the evaporation rate a„ is 
much larger and this larger evaporation impedes the for- 
mation of EHD. For the same reason the formation of 
EHD slows down. The time when EHD density reaches 
its maximum is about 10^ ps compared to about 40 ps 
under the same excitation at 12 K . The relative mean 
square of EHD distribution (An)^, which describes the 
fluctuation of droplet distribution, is very small when the 
system is in the steady state. 
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FIG. 5: Time dependence of the "exciton density" n^, the 
moments of EHD distribution, a;o,a;i,a;2, for diamond at 100 
K, the average number of e-h pairs per d rop, n and the relative 
mean square of the droplet distribution (An)^ under different 
excitation. We use Gaussian pulse excitation with = 0.1 
ps and to = 0.1 ps: (a): Go = 6 x 10^° cm'^s"^; (b): Go = 



The kinetics at a higher excitation intensity of 6.1 
mJ/cm^ is plotted in Fig. 5(b). Note that the curve of 
xi in Fig. 5(b) is comparable with that in Fig. 4(b) of 
Ref. 5 — where the excitation density and temperature are 
similar to those in the calculation except that our re- 
sult is plotted in a logarithmic scale and in a larger time 
range. Moreover the kinetics shows two main differences 
from the low excitation case: First, the average drop size 
n is much smaller: about 10^; Second, the density of all 
e-h pairs which are condensed into droplets, xi, is much 
larger around 10^* cm~^ which is at least four orders 
of magnitude larger. 

The high excitation creates a large number of excitons, 
which then produce a large number of small droplets. 
However this process, together with the growth of the 
newly created small droplets, tend to reduce n^. As 
rix becomes smaller, the thermal potential in Eq. (Al) 
forms a higher barrier between excitons/multiexciton 
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complexes and EHDs.^"^ When the thermal potential 
grows high enough, it blocks the multiexciton complexes 
to grow into small droplets and the reverse process. Thus 
for a relatively long time, the droplet density xq becomes 
stable. It is noted that xq always becomes stable before 
other moments, xi and X2, get stable as shown in Figs. 
5(a) and (b).^^'^^ Thus the formation process is separated 
into two stages: The first is the process of the growth of 
the number of droplets; and the second is the process of 
the growth of the size of the droplets. From Fig. 5(a) 
one can see that the first stage ends at about 8 x 10^ ps, 
and the second one ends at 2 x 10^ ps. The formation 
process after the first stage is a key one for the growth 
of the size of droplets. In Fig. 5(a) one can see that Ux 
stays nearly unchanged during this process, because the 
number of droplets is very small and the growth of these 
droplets requires very few excitons. In this case the ther- 
mal potential in Eq. (Al) remains nearly the same so 
that there is enough time for the e-h system to evolve 
slowly into its equilibrium, where the average drop size 
is very large, while in Fig. 5(b) one finds the process in 
80-200 ps causes to decrease by one order of magni- 
tude. The depletion of excitons prevents droplets to grow 
larger, i.e., the shortage of excitons stops the growth of 
the size of droplets when the gain £f„ and loss /„ rates 
are equal. Therefore the system can only reach a steady 
state which is in fact far away from equilibrium. It is 
noted from Eq. (7) that Uc oc (Inn^;)"^, i.e., Uc decreases 
with the increase of excitation. This, together with the 
fact that the formation rates of multiexciton complexes 
are proportional to n^, tend to increase the number of 
small droplets greatly with the increase of excitation den- 
sity. In short, a large excitation tends to create too many 
small droplets which are unable to grow into large ones 
due to the limited number of excitons. And as a result, 
the density of exciton, n^, is small while the density of 
all e-h pairs condensed in droplets, xi, is large. These 
metastable features are similar to those in Ge at a high 
enough temperature.^^ 

In summary we have studied the kinetics of EHD for- 
mation and decay at low (12 K) and high (100 K) tem- 
peratures under both low and high excitations by master 
equations. At low temperature our results are compa- 
rable with measurements reported by Nagai et al.^ The 
time evolution of exciton and EHD basically represents 
the time-resolved photo- luminescence measurement. The 
possible causes of the discrepancies between the calcula- 
tion and the experiments are discussed in detail. The 
ANPC under both low and high excitation are too small 
to form macroscopic EHDs. The smallness of ANPC and 
the time during which the system is in equilibrium indi- 
cate that the phase transition is a second order process 
as in direct semiconductors. Despite much simplification 
in the model of the master equation theory, our results 
are in good qualitative agreement with experimental re- 
sults. Our study of EHD at high temperature predicts 



that the average drop size is as large as 10®. Moreover, 
under low excitation the average size of EHDs is very 
large but the EHD density is very low. Under high exci- 
tation the average size of EHDs is much smaller, but the 
EHD density is much larger. The physics behind these 
predictions is discussed in detail. These effects demon- 
strate the metastable feature of the kinetics at high tem- 
perature in diamond. Experiments are needed to verify 
these predictions. 
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APPENDIX A: THE EXPRESSIONS OF J„<, AND 
/(net) 

The probability current from droplets in Eq. (8) is 
given by Staehli:^^ 

J) Tl J) K-L 

(Al) 

with Fx{n) — nxn^'"^ ex-'p[—^{n)/kT] denoting the dis- 
tribution function of excitons and multiexciton com- 
plexes, p = ^ 2TvkT / [d'^'^ {n) / d'n?)n^ is the width of 
the potential barrier between excitons/multiexciton com- 
plexes and EHDs. n = xi/xo is the average drop size. 

n 

C = xo/[2 J (n/n)^/^ exp(-^^^^^^^^^^) dn] is a normaliza- 

tion factor of the distribution function of EHDs and the 
thermal potential ^'(n) is given by^^ 



*(n) = -kTnln — +4TrRlan'^/^ + kTj2H^- 



a-jTd 



-)• 



(A2) 

As for /(nc, t), when > it takes the form of distri- 
bution function of excitons and multiexciton complexes: 



/(net) «F,(ne)/2 . 



(A3) 



Otherwise it takes the form of distribution function of 
EHDs: 

/(n.,i).C7(|)3/^exp(^:M^). (A4) 
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